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This  reports  the  results  obtained  from  a multiple  regression 
analysis  that  was  conducted  on  seismic  data  collected  during  field 
measurement  of  motion  originating  from  an  impulse  source.  Equations 
that  predict  the  maximum  distance  at  which  a ground  seismic  intrusion 
detector  (GS ID)  can  detect  a man  walking  are  given.  The  independent 
variable  in  these  equations  is  the  thickness  of  the  soil  overlying  a 
hard  layer. 


PRECEDING  PACK  BLANK- not 


E?r/IBONMEriTAJ.  EFFECTS  OK  MICROS  5 1 SMI  C 
WAVE  PROPA'-iATT  ON 


Background 


1.  The  term  "microseismic  wave  propagation"  as  used  herein  is  de- 
fined as  the  movement  of  waves,  whose  amplitude  is  considerably  smaller 
than  that  produced  by  an  earthquake,  through  the  earth's  crust.  More 
insight  and  knowledge  of  this  movement  may  lead  to  better  and  more  ad- 
vanced techniques  for  use  in  intrusion  detection,  security,  and  associ- 
ated areas. 

2.  Present  knowledge  has  established  the  fact  that  wave  propaga- 
tion through  the  earth's  crust  is  a very  complicated  phenomenon.  The 
earth's  crust  is,  in  general,  a heterogeneous,  dispersive,  anisotropic 
medium;  consequently,  an  exact  analytical  solution  using  the  classical 
theory  of  elasticity  for  wave  propagation  is  highly  unrealistic.  Much 
has  been  published  in  the  literature  on  wave  propagation  through  a 
homogeneous,  isotropic,  linearly  elastic  medium.  For  certain  problems, 
this  description  of  the  medium  is  adequate  and  its  use  will  yield  gener- 
alized solutions;  however,  these  solutions  are  functions  of  the  elastic 
properties  of  the  material  and  do  not  account  for  substrata  discontinu- 
ities often  encountered  within  the  earth's  crust.  Therefore,  the  solu- 
tion obtained  is  no  better  than  the  accuracy  of  the  description  of  the 
material  and  environmental  properties  that  are  used  as  input.  The  mate- 
rial must  be  described  in  terms  of  its  elastic  properties,  i.e.  the 
shear  modulus  G and  Young’s  modulus  E . Both  E and  G have  been 
expressed  as  functions  of  wave  velocity,  density,  void  ratio,  effective 
stress,  and  plasticity  index. 

3-  The  solutions  obtained  using  the  approach  mentioned  above  do 
not  adequately  describe  field  propagation  phenomena.  Consequently, 
there  is  a need  for  another  approach  to  realistically  describe  the  pa- 
rameters influencing  seismic  wave  propagation  through  near-surface  mate- 
rials. Numerous  field  experiments  using  steady-state  vibratory  and 


various  <ypes  ol'  impulse  sources  have  demonstrated  that  an  empirical  re- 
lationship could  be  developed  to  predict  the  attenuation  of  microseismic 
waves  in  layered  soil  material.  It  is  anticipated  that  this  relationship 
wiLl  be  a function  of  easily  determinable  soil  parameters. 

4.  This  study  is  a first  attempt  to  relate  environment  properties 
such  as  (a)  soil  depth,  (b)  soil  type,  and  (c)  vegetation  to  wave  propa- 
gation potential.  In  this  report  wave  propagation  potential  is  denoted 
by  the  distance  from  a given  source  at  which  a wave  can  be  detected 
using  a given  seismic  sensor.  From  the  results  obtained  it  is  hoped 
that  a more  accurate  estimate  of  a site’s  wave  propagation  character- 
istics can  be  made. 


Purpose 

5.  The  purpose  of  this  study  was  to  (a)  develop  a method  of 
analysis  that  would  yield  a model  capable  of  predicting  the  seismic  re- 
sponse of  a particular  site,  and  (b)  determine  the  existence  of  any  soil 
parameters  other  than  the  elastic  constants  which  could  be  used  to  pre- 
dict soil  response  and  seismic  wave  propagation  characteristics  at  a 
given  site.  These  parameters  were  to  be  used  as  independent  variables 
and  correlated  with  seismic  response  (maximum  detection  distance)  to 
form  a prediction  model.  This  model  was  to  contain  only  the  more  sig- 
nificant parameters  and  was  to  be  fabricated  in  a simple  and  usable  form. 

Scope 

6.  Several  techniques  were  considered  for  obtaining  field  data 
that  could  be  used  to  study  the  effect  of  soil  parameters  and  environ- 
mental conditions  such  as  temperature,  vegetation,  ground  slope,  etc., 
on  wave  propagation  through  the  soil.  It  was  decided  that  the  scope  of 
this  investigation  would  be  limited  to  seismic  field  data  originated 
from  impulse  rather  than  steady  static  sources.  This  decision  was  based 
largely  on  emphasis  the  Department  of  the  Army  has  placed  on  its  seismic 
sensor  program.  The  information  obtained  from  a study  of  this  type  has 


immediate  and  direct  application  to  the  current  work  being  done  in  the 
seismic  sensor  intrusion  detection  area. 

7.  A review  of  the  literature  produced  a minimum  of  good  field 
data  in  this  area.  This  information  was  furnished  to  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES'  by  the  Sandia  Coporation, 
Albuquerque,  N.  Mex.  These  data,  which  were  obtained  with  a ground 
seismic  intrusion  detector  (GSID)  developed  by  Sandia  Corporation,  in- 
cluded detection  distance  (for  a man  walking),  depth  of  soil,  a descrip- 
tion of  the  soil,  a description  of  the  vegetation,  site  location,  and, 
in  some  cases,  the  frequency  of  the  seismic  signal.  These  data  were 
collected  at  various  locations,  both  in  the  continental  United  States 
and  abroad.  A regression  analysis  was  conducted  on  these  data  using  a 
statistical  program  developed  by  Mr.  J.  H.  Goodnight  of  the  Experimental 
Statistics  Department  at  North  Carolina  State  University,  Raleigh,  N.  C. 
This  program  was  run  on  a GE-600  computer. 

Description  of  the  Computer  Program 

8.  The  computer  program  can  be  subdivided  into  four  parts.  The 
first  part,  part  I,  is  used  to  input  both  dependent  variables  (variables 
to  be  predicted)  and  independent  variables  (variables  used  to  predict). 
These  variables  comprise  the  measured  data  or  data  bank.  New  variables 
can  be  generated  as  functions  of  measured  variables.  For  example,  fre- 
quency and  soil  depth  might  be  measured  variables,  and  a new  variable 
(frequency  times  the  square  root  of  soil  depth)  might  be  generated.  The 
program  is  presently  limited  to  a total  of  75  variables  including  de- 
pendent and  independent  variables  (both  measured  and  generated) . 

9.  In  part  II  the  simple  statistics  of  each  variable  are  computed. 
In  this  case  simple  statistics  refers  to  the  sum,  mean,  sum  of  the 
squares,  variance,  and  standard  deviation.  Bivariant  statistics  can 
also  be  computed.  Here  each  variable  is  correlated  with  every  other 
variable  (one  at  a time).  This  analysis  helps  determine  if  measured 
variables  are  interrelated. 

10.  In  part  III  the  dependent  variable  is  specified  and  the 
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computer  searches  the  independent  variables  (first  one  at  a time,  then 
two  at  a time,  then  three,  etc.)  to  see  what  variable  and/or  combination 
of  variables  correlates  best  with  the  dependent  variable. 

11.  Based  on  the  results  of  part  III  the  dependent  variable  and 
independent  variables  are  specified;  the  computer  then  uses  the  Doolittle 
matrix  inversion  technique  to  generate  the  regression  equation.  This 
equation  is  of  the  form: 

Y = Bo  + Eixi  + V2  » •••  ’Vi 


where 


xL,x2, 

*B  ,B/}  j 


Y = dependent  variable 
.X^  = independent  variables 
= regression  coefficients 


This  equation  is  the  equation  of  best  fit  through  the  data.  "Goodness 

of  fit"  is  measured  by  the  value  of  the  squared  correlation  coefficient 
2 2 

(r  ).  If  the  r value  is  1,  perfect  correlation  has  been  obtained. 

2 

If  the  r value  is  0,  no  correlation  has  been  obtained. 


Site  Selection  and  Layout 


12.  The  data  compiled  for  use 
the  following  areas : 

Location 

West  Germany 
Fort  Hood,  Tex. 

South  Vietnam 
Baraboo,  Wis. 

Chula  Vista,  Calif. 

Panama 

Great  Britain 


in  this  study  were  gathered  from 


No.  of 

Date  of 

Sites 

Testing 

30 

Sep-Oct  70 

28 

Apr  70 

6 

Sep  70 

4 

Dec  70 

3 

May  70 

2 

Mar  70 

2 

Oct  70 

(Continued) 


Location 

Fort  Bliss,  Tex. 

Camp  Pendleton,  Calif. 

29  Palms,  Calif. 

Tonopah,  Nev. 

Eglin  Air  Force  Base,  Fla. 


No . of 
Si  tes 


Date  of 
Testing 

Nov  67 
Feb  70 
Feb  70 
Jan  70 
Oct  68 


Total  80 


The  site  locations  were  chosen  so  as  to  give  a wide  range  of  environ- 
mental characteristics.  The  test  sites  were  laid  out  in  the  following 
manner.  First,  a 120- m straight  line  was  staked  off.  This  line  was  then 
subdivided  into  twenty-four  5-m  segments  and  marked,  starting  with  zero 
at  the  center  of  the  line  and  increasing  in  both  directions  (see  fig.  1). 

60  55  50  4 5 W 35  30  25  20  15  10  5 0 5 10  15  20  25  30  35  40  45  50  55  60 


''Point  A 

(Approximately  3 * from  point  0) 


Scale:  Meters 


Fig.  1.  Typical  site  layout 


At  the  point  marked  zero,  a second  line  was  staked  off  perpendicular  to 
the  first  line.  A distance  from  the  first  line  of  approximately  3 m on 
the  second  line  was  marked.  This  location  was  identified  as  point  A, 
thus  completing  the  site  layout. 

Test  Procedure 


13.  With  the  site  layout  complete,  a GSID  was  placed  on  the 
ground  surface  at  point  A.  The  GSID  was  activated,  and  the  gain  control 
was  set  at  a specific  range  setting.  The  sensitivity  of  the  instrument 
can  be  varied  by  adjustment  of  the  gain  control.  A man  then  walked  down 
the  line,  and  the  distance  at  which  the  detector  first  indicated  ground 
motion  was  recorded.  He  continued  to  walk,  passing  station  zero,  until 
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no  motion  was  registered.  This  distance  was  also  recorded.  The  gain 
control  was  then  changed  to  other  range  settings  and  the  man  repeated 
the  walk.  In  this  manner  the  extreme  distance  at  which  a man  walking 
could  be  detected  by  a GSID  and  the  particular  gain  setting  used  were 
recorded.  In  some  cases,  additional  electronic  equipment  was  used  that 
enabled  operators  to  record  the  ground-motion  signal  and,  in  turn,  de- 
termine the  frequency  of  the  wave  being  received;  however,  this  was  done 
at  only  approximately  50  percent  of  the  sites.  The  site  environmental 
data  were  then  recorded.  These  data  included  a qualitative  description 
of  the  soil  and  of  the  vegetation.  In  addition,  the  depth  of  soil  to  a 
hard  layer  (rock  or  heavily  overconsolidated  soil)  was  determined.  If 
such  a layer  was  deeper  than  3 m,  it  was  recorded  as  "+3  m."  These 
depths  were  determined  by  both  hand  auger  and  visual  estimation.  The 
visual  estimate  was  based  on  existing  terrain  features,  air  photos,  etc. 

Presentation  of  GSID  and  Site  Data 

14.  The  GSID  and  site  data  are  presented  in  table  1.  It  should 
be  noted  that  depths  to  a hard  layer  ranged  from  0 to  +300  cm;  i.e.  any 
depth  to  rock  greater  than  300  cm  was  listed  as  +300. 

Analysis  and  Results 

15.  Sensor  response  was  emphasized  during  the  testing  program. 
Minimum  effort  was  expended  to  describe  soil  (surface  and  subsurface) 
and  vegetation  conditions.  With  this  qualitative  data,  attempts  were 
made  to  quantitatively  describe  each  site.  In  order  to  do  this  a number 
was  assigned  to  each  soil  type,  i.e.  rock  = 0,  sand  = 1.0,  silt  = 2.0, 
and  clay  = 3.0.  With  this  system  a silty  clay  would  be  a 2.5.  Soil 
homogeneity  was  also  ranked;  i.e.,  good  = 1.0,  medium  = 2.0,  and 

poor  = 3-0.  With  this  system  a uniform  sand  would  have  a homogeneity 
of  1.0  and  rocky  clay  might  be  3*0.  Vegetation  was  also  ranked  in  the 
following  manner:  no  vegetation  = 1.0,  grass  = 2.0,  trees  = 3-0,  and 

jungle  = 4.0. 
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16.  The  system  mentioned  above  was  used  to  describe  data  col- 
lected at  sites  in  West  Germany;  at  Fort  Hood,  Tex.;  Baraboo,  Wis.; 

Chula  Vista,  Calif.;  and  in  England  and  these  descriptions  were  fed  into 
the  computer.  The  variables  were  subdivided  as  follows: 

Dependent  variable 

X^=  detection  distance  of  GSID  at  gain  5,  m 


Independent  variables 

X0  = soil  depth,  m 
X^  = soil  homogeneity 
X^  = soil  type 
Xj-  = vegetation 
Xg  = frequency,  Hz 
X^  = logarithm  of  X? 


It  should  be  noted  that  when  frequency  was  not  determined  it  was  re- 
corded as  zero.  It  was  found  that  the  best  combination  of  three  vari- 
ables was  X0  , X , and  X,  . These  variables  fit  the  data  with  an 

2 ^3  't 

r value  of  0.70.  If  only  one  independent  variable  is  used,  fits 

2 

the  data  best  with  an  r value  of  0.66.  Several  other  attempts  were 
made  to  correlate  the  data  using  this  crude  quantitative  environmental 

site  description.  Since  the  addition  of  X_  and  X^  only  increases  the 

P ^ 

r value  0.04,  it  was  concluded  that  this  type  of  site  description  does 

not  significantly  improve  the  ability  to  predict  detection  distance,  and 
that  Xp  was  the  most  significant  of  the  measured  variables. 

17.  Next  the  data  were  subdivided  into  two  groups:  (a)  data 

points  with  frequency  measurements  and  (b)  data  points  without  frequency 
measurements.  The  detection  distance  obtained  at  gain  5 was  the  only  de- 
tection distance  used  in  this  analysis,  because  more  data  were  available 
for  gain  5 than  either  of  the  other  gains.  Often  no  detection  distance 
was  determinable  at  gains  4 and  3 or  no  data  were  taken  at  these  gains . 
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18.  A regression  analysis  was  conducted  on  the  data  in  group  a. 
Various  functions  of  depth  and  frequency  were  tried  as  the  independent 
variable , keeping  the  detection  distance  as  the  dependent  variable . 

From  this  analysis,  it  appears  that  the  simplest  equation  of  best  fit  is 


where 


depth  to  a hard  layer  (s  3 m) , m 


Equation  1 fits  the  data  with  a correlation  coefficient  r of  0.7115- 
This  model  is  limited  by  the  restricted  amount  of  available  data,  use  of 
only  one  independent  variable  , and  because  the  phenomenon  of  wave 
propagation  is  obviously  a function  of  other  environmental  characteris- 
tics. In  this  light,  an  r value  of  0.7115  is  considered  both  reason- 
able and  acceptable. 

19-  Table  2 shows  the  detection  distances  predicted  by  equation  1, 
the  actual  observations,  and  the  differences  between  the  predicted  and 
observed  distances. 

20.  If  equation  1 is  divided  by  two  and  rounded  off  to  the  nearest 
whole  number , i . e . , 


it  would  yield  a conservative  value  of  approximately  85  percent  of 

the  time.  By  conservative,  it  is  meant  that  the  predicted  detection 
distance  is  less  than  the  observed  detection  distance.  Equation  la  is 
not  conservative  for  sites  10,  35,  60,  66,  and  74  (see  table  2).  Next, 
as  a check,  equation  1 was  used  to  predict  values  from  the  data  in 
group  b,  which  were  not  used  to  generate  this  equation.  Table  3 is  a 
tabulation  of  typical  results  of  this  exercise.  The  analysis  indicates 
that  equation  1 yields  reasonable  prediction  about  60  percent  of  the 
time.  However,  if  the  predicted  values  are  divided  by  2 they  are  less 
than  the  observed  value  (conservative)  most  of  the  time.  Therefore 
equation  la  can  be  used  to  conservatively  predict  the  distance  at  which 


the  GSID  can  detect  a man  walking.  Since  Xp  can  never  be  negative, 
the  minimum  detection  distance  pz-edicted  by  equations  1 and  la  is  7*57 
and  4 meters,  respectively. 

21.  A second  regression  analysis  was  conducted,  this  time  using 
all  the  data  (groups  a and  b).  This  analysis  showed  that 

XL  = 7.38  + 3-97  (2) 

is  the  simplest  equation  of  best  fit.  This  equation  fits  the  data  with 
a correlation  coefficient  of  0.772.  Table  4 presents  the  predicted 
values  using  equation  2,  the  observed  values,  and  differences.  A study 
of  table  4 indicates  that  equation  2 predicts  detection  distances  within 
25  percent  of  observed  detection  distances  approximately  75  percent  of 
the  time,  if  sites  with  observed  detection  distances  of  less  than  6 m 
are  neglected. 

22.  If  equation  2 is  reduced  by  one-half,  then 

Xx  = | (7.38  + 3-97  Xp)«  4.0  + 2.0  Xp  (2a) 

If  this  modification  is  applied  to  table  4,  equation  2a  gives  a con- 
servative prediction  (less  than  observed)  approximately  85  percent  of 
the  time. 

23.  All  the  equations  presented  in  this  report  are  parabolas. 

This  indicates  that  the  detection  distance  varies  as  the  square  of  the 
soil  depth.  When  the  thickness  of  soil  is  small  (approaches  0) , equa- 
tions la  and  2a  predict  a constant  detection  distance  of  approximately 

4 m.  Since  the  soil  thickness  can  never  be  negative  these  equations  are 
unable  to  predict  detection  distances  less  than  4 m. 

Conclusions 


24.  It  should  be  emphasized  that  the  conclusions  drawn  in  this 
report  are  based  on  limited  data.  The  following  are  the  primary  con- 
clusions, as  indicated  by  the  results  of  this  study. 

a.  If  the  subsurface  conditions  (depth  to  rock)  to  a depth 
of  3 m are  known,  a conservative  estimate  of  the  detec- 
tion distance  can  generally  be  made.  This  detection 
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distance  is  a measure  of  the  site's  ability  to  propagate 
a signal. 

b.  A mathematical  model  developed  to  predict  the  distance  at 
which  a GSID  seismic  sensor  will  detect  a man  walking  is 
least  reliable  when  the  soil  layer  thickness  Xp  is 
small.  The  minimum  predictable  detection  distance  exists 
when  Xp  = 0.  For  this  case,  Xq  « 7*5  if  regression 
equations  1 and  2 are  being  used  and  Xq  « 3-75  if  equa- 
tions la  and  2a  are  used. 

c.  It  appears  that  the  distance  at  which  a signal  can  be  re- 
ceived by  a seismic  geophone  varies  as  the  square  of  the 
depth  of  soil  above  a hard  layer. 


Recommendation  for  Future  Research 


25.  As  a result  of  this  investigation,  it  is  believed  that  if 
data  were  available  which  yielded  more  information  about  the  site 
environment,  a better  and  more  reliable  model  could  be  developed.  It  is 
recommended  that  data  be  collected  in  which  the  following  are  measured: 

a.  Water  content 

b.  Density 

£.  Void  ratio 

d.  Cone  index 

e.  Layer  thickness 

f . Compression  and  shear  wave  velocities 

£.  Predominant  frequency 

h.  Background  noise 

i.  Temperature  (soil  and  air) 

Wind  speed  and  direction 

k.  Vegetation  conditions 

l.  Detection  distance — from  a footstep  and  a controlled 
energy  source 

The  seismic  information  should  be  obtained  with  calibrated  instrumentation 
and  trained  personnel.  Presently,  the  WES  is  obtaining  data  that  will 
include  these  parameters.  It  is  recommended  that  a similar  study  be 
conducted  incorporating  these  parameters  in  a multiple  regression 
analysis. 
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Appropriation  No.  21X2OU0  108-3508  P50IB 


'This  reports  the  results  c:  ^nined  frm  a multiple  regression  analysis  that  was  con- 
ducted on  seismic  data  collec*  ;d  during  field  measurement  of  motion  originating  from 
an  impulse  source.  Equations  that  predict  the  maximum  distance  at  which  a ground 
seismic  intrusion  detector  (GSID)  can  detect  a man  walking  are  given.  The  indepen- 
dent variable  in  these  equations  is  the  thickness  of  the  soil  overlying  a hard  layer 


irity  Classification 


DOCUMENT  CONTROL  DATA  R & D 

Security  ( la*  ai  He  at  ton  ot  till*,  body  ol  abstract  and  mdamlng  annotation  muat  ba  antarad  whan  tha  ovarall  raport  /«  clmaalllmdj 

I ORI6INATINQ  ACTIVITY  ( Corporals  author) 

la.  REPORT  EBCURlTT  CLAMlPIC  ATION 

U.  S.  Army  Engineer  Waterways  Experiment  Station 

Unclassified 

Vicksburg,  Miss. 

26  CROUP 

» REPORT  TlTLt 

^ENVIRONMENTAL  EFFECTS  ON  MICROSEISKIC  WAVE ^PROPAGATION 

IT 

10  Ol  ITRl  BUTlON  STATEMENT 

Approved  for  public  release;  distribution  i 

unlimited. 

ii  Supplement a«> 

r notes 

12.  SPONSORING  MILITARV  ACTIVITY 

Office,  Chief  of  Engineers,  U.  S.  Army 

| ■■ 

Washington,  D.  C. 

Unclassif ied 

Security  Classification 


Environmental  effects 
Microseismic  waves 
Regression  analysis 
Seismic  detection 
Wave  propagation 


Unclassified 


